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ABSTRACT 

The  multicomponent  diffusion  coefficient  Du  „  has  been  evaluated  from 

H,h2 

measurements  of  the  relative  atom  density  in  a  semi-infinite  catalytic 
cylinder  as  a  function  of  (a)  distance  from' the  atom  source  under  steady- 
state  conditions,  and  ( b)  time  during  the  transient  condition  associated 
with  removal  of  the  atom  source  after  attainment  of  steady  state.  An 
electron'  spin  resonance  spectrometer  was  employed  for  -atom-density 
measurements  and  a  radio- frequency  discharge  for  atom  production  in  the 
total  gas-pressure  range  from  0.85  to  2.6  torr.  The  value  for  the  dif¬ 
fusion  coefficient  was  found  to  be  D  _  =  1.8  ±  0.13  cm2*s«c-1  at  1 

H,  n2 

atmosphere  and  293  ±  3°K.  Experimental  data  for  the  variation  of  D„ 

«,  h2 

as  a  function  of  temperature  (293  to  719°K)  are  analyzed  in  terms  of 
different  interaction  potentials  for  the  atom-molecule  collision. 


*Thi3  work  was  sponsored  by  Project  Squid  which  is  supported  by  the 
office  of  Naval  Research,  Department  of  the  Navy,  under  Contract 
Nonr  3623(00),  NR-098-038. 


INTRODUCTION 


A  number  of  publications1-?  2> 3  have  dealt  with  the  theoretical  calcula¬ 
tion  of  the  binary  diffusion  coefficient  of  atomic  hydrogen  in  molecular 
hydrogen  (D  ).  Also  an  estimate  of  this  parameter  has  been  made  from 
viscosity  measui;ements4-> 5  on  the  basis  of  theoretical  interaction 
potentials.  However  no  direct  determination  of  the  binary  diffusion 
coefficient  has  been  reported  to  ascertain  the  degree  of  validity  of  the 
theoretical  analyses,  which  involve  calculations  based  on  multiple  inter¬ 
action  potentials  for  species  with  unfilled  electronic  shells.* 

In  a  preceding  publication  we  demonstrated  the  feasibility  of  measuring 

diffusive  transport  coefficients  by  a  combination  of  steady-state  and 

transient  experiments.7  We  have  applied  the  principle  of  this  method  to 

the  study  of  D  . 

h,h2 

APPARATUS 

The  diffusion  coefficient  of  hydrogen  atoms  in  hydrogen  was  measured 
by  means  of  the  apparatus  shown  schematically  in  Fig.  1.  The  central 
part  of  the  vacuum  apparatus  is  a  quartz  diffusion  tube  9.4  mm  i.d.  which 
includes  a  constriction  2  mm  i.d.  and  3  mm  long.  On  one  side  of  the  con¬ 
striction  the  tube  (45  cm  in  length)  passes  through  the  ESR  cavity.  The 
quartz  section  is  attached  by  means  of  an  "0"  ring  gland  to  a  pyrex  tube 
(10  mm  i.d.  and  50  cm  in  length)  on  the  end  of, -which  is  a  pressure  gauge 
.(Autovac  pirani  tube).  On  the  other  side  of  ..'he  constriction  the  quartz 
-tube  passes  through  an  rf  discharge  coil  and  a  liquid  nitrogen  trap 
(6  mm  i.d.,  50  cm  from  the  constriction).  This  end  of  the  apparatus 
communicates  with  the  rest  of  the  vacuum  system  comprising  metal  valves, 
a  McLeod  gauge  for  pressure  calibration,  and  the  supply  of  hydrogen. 

The  hydrogen  (Matheson  &  Company,  prepurified  grade,  purity  >99.95$) 
was  passed  through  a  de-oxo  purifier  and  two  traps  cooled  by  liquid 
nitrogen  before  entering  the  apparatus.  For  cleaning, .the  glass:. tubing 
was  exposed  successively  to  concentrated  nitric  acid,  concentrated 
hydrochloric  acid,  10$  hydrofluoric  acid,  and  chromic  ^acid-sulfuric  acid 
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cleaning  solution.  ^Finally  the  tube  was  rinsed  with  large  quantities  of 
distilled  water.  The  last  solution  was  tried  in  an  attempt  to  change 
the  activity  of  the  quartz  walls  for  hydrogen  atom  recombination.  The 
treatment  resulted  in  a  decrease  of  catalytic  activity  by  a  factor  of 
about  two. 

The\r»f  discharge  was  .excited  by  inductive  coupling  to  an  18  mc/sec 
trarsniitter,  and  the  region  of  the  rf  coil  was  cooled  by  a,  fan.  The 
center  of  the  rf  coil,  around  the  quartz  tube,  was  located  about  7  cm 
from  the  constriction  and  about  25  cm  from  the  center  of  the  magnet. 

This  configuration  provided  satisfactory  atom  densities  in  the  diffusion 
tube.  Also  it  prevented  the  plasma  from  extending  beyond  the  constric¬ 
tion.  At  low'  pressure  and  low  rf  power,  a  Tesla  coil  was  required  to 
initiate  the  discharge  for  the  'transient  experiments. 

For  the  transient  measurements  the  pulsing  of  the  discharge,  the  initia¬ 
tion  of  the  Tesla  coil,  and  the  triggering  of  the  time-averagihg.eogiputer 
(Varian  C1024)  were  synchronized.  The  time  basis  for  these  measurements 
was  obtained  from  an  electromechanical  timer  ( FI sxopulse,  Eagle  Signal 
Corp.).,  with  which  ON-  and  OFF-tirae  of  the  discharge  could  be  adjusted 
independently.  This  timer  operated  a  mercury  relay,  one  pole  of  which 
keyed  the  transmitter.  Another  pole  of  the  mercury  relay  operated  an 
auxiliary  circuit  that  turned  the  Tesla  coil  ON  for  0.1  sec  at  the  time 
the  transmitter  was  turned  ON.  Operation  of  the  Tesla  coil  was  limited 
to  as  short  a>  period  as  possible  in  order  to  avoid  electrical  interfer¬ 
ence  with  the'sensitive* electronic  measurement  devices.  Since  the  surface- 
catalyzed  recombination  of  atomic  hydrogen  is  a  reaction  of  first-order 
in  gas  atom  density,  measurements  of  relative  atom  densities  are  adequate 
for  the  evaluation  of  the  surface  activity  and  the  diffusion  coefficient. 

For  determination  of  relative  hydrogen  atori  density  and  its  rate  of 
decay  in  the  transient  experiments,  we  employed  an  electron  spin  resonance 
(ESR)  spectrometer  (Varian  V-4502,  X-band)  equipped  with  a  12-in.  magnet 
and  a  time-averaging  computer.  The  magnetic  field  was  controlled  by  a 
Fieldial  (Varian  Mark  I)  for  the  steady-state  measurements^  Since  the 
Hall  effect  sensor  of  the  Fieldial  was  sensitive  to  the  rf  field  of  the 
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discharge,  current-regulation  of  the  magnetic  field  was  needed  during 
transient  measurements.  A  rectangular  cavity  (TE102)  was  attached  to 
the  low  power  connection  of  the  microwave  bridge.  The  microwave  power 
was  adjusted  for  maximum  signal  intensity.  A  modulation  frequency  of 
100  kc/sec  was  employed,  and  the  modulation  amplitude  was  adjusted  for 
maximum  signal  intensity  which  also  provided  substantial  line  broadening 
(8  Oe)  so  that  during  the  transient  measurement  any  fluctuations  in  the 
magnetic  field  were  not  significant. 

A  chip  of  ruby  was  affixed  to  one  inside  wall  of  the  cavity  to  provide 
a  monitor  for  its  sensitivity  as  it  was  moved  along  the  quartz  tube. 
Movement  of  the  cavity  was  achieved  by  simply  sliding  the  microwave 
bridge,  with  its  attached  waveguide  leading  to  the  cavity,,  on  the  hori¬ 
zontal  table  mounted  above  the  magnet.  This  procedure  permitted  a 
horizontal  cavity  movement  of  7  cm  with  a  variation  of  reference  signal 
Intensity  of  less  than  1.5$;  greater  movement  resulted  in  line  distortion* 
due  to  Inhomogeneous  magnetic  fields. 

EXPERIMENTAL  DETAILS  AND  RESULTS 
A.  Transient  Measurement 

The  transient  experiments  Involved  study  of  the  rate  of  atom  density 
decay  from  the  steady-state  value  at  a  given  point  in  the  cylindrical 
reactor  after  rapid  extinction  of  the  atom  source For  the  determination 
of  the  atom  density,  the  magnetic  field  was  adjusted  for  a  peak  of  the 
ESR  derivative  signal  and  the  transient  decay  was  time-averaged  by  cycl¬ 
ing  the  discharge  with  an  ON- time  of  2  sec  and  an  OFF-time  of  1  sec. 
Collection  of  the  experimental  data  was  begun  after  exposure  of  the 
apparatus  to  the  discharge  cycling  for  a  period  of  15  min.  At  a  given 

pressure  the  power  input  into  thf /  rf  discharge  was  adjusted  to  a  level 

adequate  for  atom  detection,  yet  low  enough  to  make  the  atomic  species 
the  minor  component  of  the  gas  mixture.  The  cavity  was  situated  as  close 
to  the  constriction  as  possible  in  order  to  obtain  the  highest  atom 
density.  However,  the  distance  between  the  cavity  center  and  the  con¬ 
striction  could  not  be  less  than  15  cm,  a  distance  dictated  by  a  compromise 
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between  the  interaction  of  the  magnet  with  the  plasma  and  the  limitations 
of  magnetic  field  homogeneity.  In  order  to  obtain  an  adequate  signal-to- 
nolse  ratio  for  the  transient  measurement  of  hydrogen  atom  decay,  the 
results  of  about  100  experiments  were  averaged  on  the  computer.  The 
electrical  time  constant  associated  with  the  ESR  signal  to  the  computer 
was  adjusted  to  1  msec,  which  is  less  than  5$  of  the  time  for  l/e  decay 
of  the  signal;  the  decay  rate  was  unaffected  by  shorter  electrical  time 
constants. 

Typical  transient  decays  of  logarithm  of  the  relative  atom  density 
versus  time  are  shown  in  Fig.  2  for  several  pressures.  For  convenience 
the  atom  densities  in  the  figure  have  been  normalized  at  time  zero 
because  the  relative  atom  densities  depended  not  only  on  the  pressure 
but  also  on  the  rf  power,  which  was  adjusted  to  a  convenient  value. 

B.  Steady-State  Measurement 

For  the  determination  of  the  atom  density  gradient  along  the  length  of 
the  quartz  diffusion  tube,  the  rf  discharge  was  operated  continuously 
and  the  atom  density  was  recorded  for  various  positions  of  the  cavity 
along  the  tube.  At  a  given  pressure, measurements  were  made  after  the 
atom  density  at  some  position,  x,  was  constant.  Usually  a  steady-state 
condition  occurred  after  about  10  mins. 

The  steady-state  data  obtained  at  the  pressures  corresponding  to  the 
transient  measurements  of  Fig.  2  are  plotted  in  Fig.  3.  The  linear 
relationships  between  logarithm  of!  atom  density  and  distance  x  at  given 
pressures  satisfy  the  requirements  of  the  theory  for  an  infinite  tube 
with  catalytic  walls  for  which  the  heterogeneous  atom  recombination  is 
first  order  in  atom  density.  The  linearity  of  the  curves  further  estab¬ 
lishes  uniformity  of  the  wall  activity  along  the  tube.  The  magnitude  of 
the  slopes  are  labeled  on  each  curve  (Fig.  3)  .  They  are  found  to  be  pro¬ 
portional  to  the  square  root  of  the  total  pressure  in  accordance  with 
theory.8 
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THEORETICAL  ANALYSIS 


The  mathematical  analysis  requires  the  steady-state  and  transient 
solutions  for  a  reacting  system  involving  a  cylinder  of  infinite  length 
and  diffusion  of  a  reactant  (present  at  low  concentrations)  as  a  result 
of  a  catalytic  reaction  by  first-order  heterogeneous  kinetics.  The 
transport  equation  in  cylindrical  coordinates  applicable  to  this  prob¬ 
lem  reads9: 


n 

r 


+  n  ) 
xx 


U) 


where  D  is  the  diffusion  coefficient,  n  the  atom  density  (in  ihe  gas 
phase),  and  r  the  radical  coordinate  in  the  cylindrical  reactor  of 
radius  R. 

The  boundary  conditions  are: 


n(  r,  o,  t)  =o  for  t  >  o 

at 

X 

if 

o 

(2) 

n  -•  o 

as 

x  -•  00 

(3) 

ir(o,x,  t)  =  o 

at 

r  =  o 

(4) 

+  jjf-  n(R,x,t)  =  o 

at 

r  =  R, 

(5) 

where  6  represents  a  dimensionless  parameter  involving  the  catalytic 
efficiency  Y  of  the  cylinder  walls 

6  =  [4D/(YCR)]C(2-Y)/2l  ,  (6) 

V 

with  0  Y  ^  and  C,  the  mean  atomic  velocity  of  the  reactant. 


9The  subscripts  denote  derivative  with  respect  to  (he  axial  direction 
(x),  the  radial  (transverse)  direction  (r),  and  time  ( t) . 
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The  boundary  condition  given  by  Eq.  (2)  specifies  the  existence  of  a 
strong  "sink"  in  the  region  of  the  atom  source  at  t  >  o,  a  condition 
prevailing  during  our  experimental  measurements.  The  quartz  walls  warmed 
by  exposure  to  the  rf  discharge  and  the  liquid  nitrogen’,  trap  in  close 
proximity  to  the  discharge  act  as  efficient  atom  sinks14'  and  cause  rapid 
decay  in  atom  density  as  soon  as  the  discharge  is  turned  off. 

The  initial  condition  for  the  time-dependent  solution  is  given  by  the* 

steady-state  atom  distribution  in  the  cylinder8 

00 

n*(r,x)  =  2^  JjSj*)  exp(-?ix)  ,  (7) 

where 

=  2  no/{R-g*d  +  62R2?2)  J^R)}  ,  (8) 

with  nQ  the  atom  density  at  x  =  o,  and  (i  =  1,2,3.  .),  the  positive 

roots  of  the  Bessel-function  equation 

J0(5R)  =  6-  R  5  jx(5r)  .  (9) 

By  mathematical  methods  similar  to  those  described  in  reference  7  one 
finds  the  transient  solution 

00 

n(r,x,t)  =  \  ^  A1P1(x,t)  r)  ,  (10) 

with  defined  by  Eq.  (8)  and 

F.(x,.t)  =  *-exp(g  x)  erfc  C§  V5t"  + — —  1 
1  11  2-tfDt 

+  exp(-5.x)  erfc  Eg.  -/Dt - .  (ll) 

1  1  2  VSt 

The  value  of  n  averaged  over  the  cross  section  of  the  tube  is  found  to  be: 

o> 

n(x,t)/nQ  =  2  F1(x,t)/r.R2?*  (l  +  fiVg*)]  .  (12) 
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For  conditions  where  the  one-term  approximation  of  Eq.  ( 12)  is  applicable 
and  the  value  of  6  is  large  (low  catalytic  activity  and/or  high  diffusion 
coefficient)  Eq.  ( 12)  reduces  to 

n(x,t)/nQ  =  |  {-exp(-g-)2  erfc  (Vlft  +  -) 

+  exp[ -(-—-) 2  ]  erfc  (VSt - - - )}'  ,  (13) 

2VDt  J 

l 

where  the  quantity  =  (2/6R2)5  =  K/D,  with  K  representing  the  rate 
constant  for  the  first-order  surface  reaction  (K  =  2D/6R2  «=  Vc/2R) .  Also 
it  is  found  that  at  long  times,  t  »  l/K,  and  large  values  of  the  argu¬ 
ments  in  the  complementary  error  function,  the  relative  reactant  decay 
becomes  an  exponential'  function  of  time 

n(x,  t)/n  [x/(KtVrnDt)  exp  (-Kt)  .  (14) 

Similarly  for  moderately  active  cylinder  walls  the  steady-state  solution 
yields 

n*(?;)/no  -  exp  C-(K/D)a]  x  .  (15) 

DISCUSSION 

The  preceding  theoretical  analysis  shows  that  the  combination  of 
transient  and  steady-state  measurements  yields  the  parameter  D^.  In 
the  system  under  examination,  the  experimental  results  for  the  variation 
of  atom  density  with  time  (after  termination  of  reactant  input  into  the 
system)  and  with  distance  ( from  the  source)  were  analyzed  by  means  of 
Eqs.  ( 13)  and  (15).  The  condition  t  »  l/K,  which  applies  to  the  approxi¬ 
mate  solution  expressed  in  Eq.  (14),  would  have  required  atom-density 
measurements  at  t  »  1  sec  with  atom  densities  so  small  as  to  make  the 
detector  output  from  the  ESR  signal  too  low  relative  to  the  noise  level* 
of  the  apparatus. 

In  the  calculation  of  the  binary  diffusion  coefficient  we  use  the 
value  of  K/D  obtained  from  the  steady-state  measurements  and  compute  the 
atom  density  for  each  of  a  number  of  values  of  Dt,  i.e.,  n  as  a  function  of 
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Dt,  by  means  of  Eq.  (13).  The  transient  experiment  yields  the  time  of 
occurrence  of  each  of  these  atom  densities  n(t),  which  by  substitution 
into  the  n(bt).  values  provide  the  parameter  D.  The  values  of  D  so 
obtained  for  a  series  of  points  on  the  experimental  n(  t)  curve  exhibit  a 
deviation  of  less  than  10$  from  the  average  value.  A  typical  example  of 
the- calculated  transient  atom  density  employing  an  average  value  of  D  is 
shown  in  Fig.  4.  The  diffusion  coefficients  derived  from  the  measure¬ 
ments  <at  the  various  total  gas  pressures  are  summarized  in  Table  I.  The 

average  value  of  D  at  '293  ±  3°K  and  1  atmosphere  is  found  to  be 
H,Ha 

1.80  ±0.13  cm* •sec-1. 

The  relative  changes  in  Du  as  a  function  of  temperature  reported  in 

reference11  may  now  be  converted:  into  absolute  values  based  on  results 

presented  in  this  study.  The  data  span  a  temperature  range  from  293°K 

to  719°K  (Table  II).  They  allow  comparison  with  various  theoretical 

analyses  of  the  transport  coefficients  of  the  H,Ha  system.  The  data 

shown  in  Fig.  5  compare  the  experimental  diffusion  coefficients  with 

recent  semiemplrlcal  calculations  by  Weissman  and  Mason3  on  the  basis 

of  a  modified  Buckingham  interaction  potential  (exp-6;  e/k  =  16.7°K; 
o 

r  a  3.847  A;  and  or  *  12.45).  The  theoretical  values  of  D  fall  on  a 
m  h,  Hj 

curve  (dashed  line.  Fig.  5)  well  above  the  experimental  one  at  elevated 

temperatures.  By  adjusting  the  magnitude  of  the  parameter  r  to  a  value 
o  m 

of  3.90  A  in  the  (exp-6)  potential  function  (solid  line,  Fig.  5),  much 

closer  agreement  can  be  obtained.  Also  it  should  be  noted  that  the 

measured  value  of  D  at  293°K  is  in  satisfactory  agreement  with  one 
H,  H2 

derived  from  recent  viscosity  measurements9  for  the  H,H2  system  and 
analyzed  in  terms  of  the  same  (exp-6)  interaction-potential. 

The  experimental  results  obtained  under  steady-state  conditions  (Fig.  3) 
lend  themselves  to  an  evaluation  of  the  catalytic  surface  activity  of 
the  quartz  walls  for  hydrogen-atom  recombination.  The  steady-state  axial 
atom  decay  is  an  exponential  function  of  the  distance  from  the  source 
(Eq.  15)  and,  oil  a  serailogarithmic  plot,  the  slope  is  proportional  to 
(K/D)^,  An  analysis  of  the  data  is  summarized  in  Table  III.  The  value 
of  V  B  1.8  x  10-8  for  quartz  is  considerably  lower  than  that  reported 


9 


for  Pyrex, 10  but  it  is  of  the  same  magnitude  as  that  for  silica  at  room 
temperature. 13 
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Table  I.  Experimental  values  of  binary  diffusion  coefficient 
in  H,Ha  system  at  293  i  3°K. 


Total 

pressure 

torr 

Dh,h2 

cm2 ‘sec-1 

D  *P 

h,h2 

atm* cm2 •sec"1 

0.85 

1540 

1.73 

1.43 

795 

1.50 

2.10 

698 

1.92 

2.15 

590 

1,67 

2.56 

530 

1.81 

^  Av.  =  1.80  dk  0.13^ 
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Table  II.  Diffusion  coefficients  of  the  H,HS  system 
at  elevated  temperatures. 


Temp 

°K 

V 

(atm*cm9 

relative8 

•  P 

i 

•sec-1) 

absolute 

293 

1.00 

1.80b 

349 

1.29 

2.32 

372 

1.67 

3.01 

375 

1.42 

2.56 

422 

1.76 

3.16 

455 

2.25 

4.05 

498 

2.14 

3.86 

518 

2.60 

4.68 

666 

3.10 

5.58 

719 

4.10 

7.39 

aH.  Wise,  J.  Chem.  Phys.,  34,  2139  (1961). 
bthis  paper. 
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Table  III.  Catalytic  surface  activity  of  quartz  for  hydrogen  atoms 

at  293°K. 


’  p 

( torr) 

D 

(cma*sec-1) 

l 

(K/D)5 

(cm) 

K 

(sec-1) 

Y(x  108) 

0.85 

1540 

0.060 

5.5 

2.0 

1.43 

7S5 

0.076 

4.6 

1.8 

2.10 

698 

0.086 

5.0 

1.9 

2.15 

590 

0.085 

4.3 

1.6 

2.56 

530 

0.092 

4.1 

1.6 

(  Av,  e‘1.8  / 
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